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^ , Abstract 

A systematic way of construction of (2+l)-dimensional dispersionless inte- 
grable Hamiltonian systems is presented. The method is based on the so-called 
central extension procedure and classical R-matrix applied to the Poisson algebras 
of formal Laurent series. Results are illustrated with the known and new (2+1)- 
' dimensional dispersionless systems. 

(To appear in J. Phys. A: Math. Gen.) 

1 Introduction 

O 

Dispersionless integrable Hamiltonian systems are often considered as a quasi-classical 
limit of the related soliton systems (Takasaki and Takebe JIJ , Konopelchenko and Alonso 
and the literature quoted there). Nevertheless, it seems that a more systematic 
approach, allowing a construction of such systems from scratch, is necessary. Actually, 
we are interested in a systematic way of construction of a class of dispersionless systems 
having a Hamiltonian structure, and infinite hierarchy of symmetries and conservation 
laws. One method of doing it is based on the classical i?-matrix theory. As well known, 
the i?-matrix formalism proved very fruitful in a systematic construction of soliton 
systems (see for example ||-|| and the literature quoted there). So, it seems reasonable 
to develop such a formalism for dispersionless systems. Recently, an important progress 
in that direction was made by Luen-Chau Li 0. In paper 0] we apply his results 
to a particular class of Poisson algebras || in order to construct multi-Hamiltonian 
(l + l)-dimensional dispersionless systems. 

Having such an effective theory for constructing multi-Hamiltonian dispersionless 
dynamical systems in (l+l)-dimensions, we were prompted to extend this method onto 
(2+l)-dimensions. The central extension was considered in early works by Reyman and 
Semenov-Tian-Shansky [[], |TD[ and also by Prykarpatsky [II], [EJ. The central extension 



approach to integrable field and lattice- field systems was presented as well in [13, [14]. 



*E-mail: blaszakm@amu.edu.pl 
^E-mail: bszablik@amu.edu.pl 



2 



Maciej Blaszak and Blazej M. Szablikowski 



As our construction leads in general to nonlocal equations, we will understand by a 
dispersionless systems in (2+l)-dimension PDEs of the form 

r\ 71 r\ 71 r\ 

3=1 3=1 U 

where Vij and Wij are pseudo-differential operators of formal symbols T> = d~ x d y . 

The paper is organized as follows. In section 2, we briefly present a number of 
basic facts and definitions of Hamiltonian dynamics on Poisson algebras concerning 
the formalism applied. In section 3, we present the general formulation of the central 
extension procedure on Poisson algebras. In section 4 and 5, we apply this and the R- 
matrix procedure to the Poisson algebras of formal Laurent series. Then in section 6 we 
illustrate our results with the known and new integrable Hamiltonian (2+l)-dimensional 
dispersionless dynamical systems. 



2 Hamiltonian dynamics on Poisson algebras: 
/^-structures 

Here, we repeat some basic facts presented in Part I to make the paper self consistent. 
The reader familiar with Part I may skip this section. 

Definition 2.1 Let A be a commutative, associative algebra with unit 1. If there is a 
Lie bracket on A such that for each element a £ A, the operator ad a : b i— > [a, b] is a 
derivation of the multiplication, then (A, [.,.]) is called a Poisson algebra. 

Thus, the Poisson algebras are Lie algebras with an additional associative algebra struc- 
ture (with commutative multiplication and unit 1) related by the derivation property 
to the Lie bracket. 

Let A be a Poisson algebra, A* the dual algebra related to A by the duality map 

A* x A -> R : (a, a) h-> {a, a), a e A, a E A* (2.1) 

and V(A*) := C°°(A*) be a space of C°°-functions on A*. Let F E V(A*), then a map 
dF : A -> A such that 



Jt F{L + w - 



= {L',dF(L)), L, V e A'', (2.2) 

t=0 

is a gradient of F. 

We confine our further considerations to such Poisson algebras A for which its dual 

A* can be identified with A. So, we assume the existence of a product (•, -)a on A which 
is symmetric, non-degenerate and ad-invariant: 

(ad a 6, c) A + (b, &d a c) A = 0, a,b,cEA. (2.3) 

Then, we can identify A* with A, (A* = A) by setting 

(a, b) = (a, b) A , a,b E A, a E A*, (2.4) 



where a E A* is identified with a E A. 
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Definition 2.2 A linear map R : A —> A is called a classical i?-matrix if the R- 

bracket 

[a,b] R :=[Ra,b] + [a,Rb}, a,b G A, (2.5) 
is a second Lie product on A. 

Theorem 2.3 j^J Let A be a Poisson algebra with Lie bracket [•, •] and non- degenerate 
ad-invariant pairing {-,-)a with respect to which the operation of multiplication is sym- 
metric, i.e. (ab,c)A = {a>, be) a, Va, b, c G A. Assume R G End(A) is a classical R- 
matrix, then for each integer n ^ — \, the formula 

{H, F} n = (L, [R(L n+1 dF), dH] + [dF, R(L n+1 dH)]) A , (2.6) 

where H, F are smooth functions on A, defines a Poisson structure on A. Moreover, all 
{■, - } n are compatible. 

The related Poisson bivectors ir n are given by the following Poisson maps 

n n :dH^ -ad L R(L n+1 dH) - L n+1 R* (ad L dH) , n ^ -1, (2.7) 

where the adjoint of R is defined by the relation 

(a,Rb) A = (R*a,b) A . (2.8) 



Notice that the bracket ( |2.6| ) with n = — 1 is just a Lie-Poisson bracket with respect 
to a Lie bracket (12. 1 



{H,F}- 1 = (L,[dF,dH) R ) A . (2.9) 

We will look for a natural set of functions in involution w.r.t. the Poisson brackets 
. A smooth function F on A is ac?-invariant if dF G ker adi, i.e 

[dF,L]=Q, LeA, (2.10) 

which are Casimir functionals of the natural Lie-Poisson bracket. 
Hence, the following Lemma is valid 

Lemma 2.4 Smooth functions on A which are ad-invariant commute in {•, •}„. The 
Hamiltonian system generated by a smooth ad-invariant function C(L) and the Poisson 
structure {•, •}„ is given by the Lax equation 

L t = [R{L n+1 dC),L], LeA. (2.11) 

For any i?-matrix each two evolution equations in the hierarchy (|2.11 ) commute due 
to the involutivity of the Casimir functions C q . Each equation admits all the Casimir 
functions as a set of conserved quantities in involution. In this sense we will regard 
Q2.11|) as a hierarchy of integrable evolution equations. 

Let us assume that an appropriate product on Poisson algebra A is given by the 
trace form tr : A —>■ R 

{a, b) A = tr (ab), a,b G A. (2.12) 
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To construct the simplest /^-structure let us assume that the Poisson algebra A can 
be split into a direct sum of Lie subalgebras A + and v4_, i.e. 

A = A+®A_, [4,4]c4 (2.13) 

Denoting the projections onto these subalgebras by P±, we define R- matrix as 

R= l -{P + -P_) (2.14) 

which is well defined. 

Two following Lemmas [n§,|| are useful in calculating of Hamiltonians H(L) from 



the gradients dH(L) 

Lemma 2.5 (Poincare). If M is a linear space, or more generally is of the star shape 
(Vl£A/{AL : ^ A ^ 1} C M), each closed k-form is exact. 

Lemma 2.6 Let M fulfill the condition of the Poincare Lemma. Then for an exact 
one- form ^{L) 

H(L) = [ ( 7 (AL),L) dX (2.15) 



is a zero-form such that dH(L) = ^{L). 

Following the above scheme, we are able to construct in a systematic way integrable 
multi-Hamiltonian dispersionless systems, with infinite hierarchy of involutive constants 
of motion and infinite hierarchy of related commuting symmetries, ones we fix a Poisson 
algebra. 

3 Central extension approach 

Assume now that the Poisson algebra A depends effectively on an independent parameter 
y £ S 1 , which naturally generates the corresponding current operator algebra C(A) = 
C°°(S l ,A) with the following modified Tr-operation: 



Tr(a) := / tr(o)dy, (3.1) 
where tr ( |2.12| ) operation is defined for the Poisson algebra A. The scalar product reads 



(a, b) c( A) ■= Tr(a • b) (3.2) 

for a and b e C(A). The current Poisson algebra C(A) can be naturally extended via the 
central extension procedure: C(A) — > C(A) = C(A) ®C with the following Lie product: 

[(a, a), (6, /3)] := ([a, b],u 2 (a, &)), (3.3) 

where a,/3 G C and uo 2 : C(A) x C(A) — > C is the standard Maurer-Cartan two-cocycle 
on C(A): 

oj 2 (a,b):=f (a, dy = Tr(a • b y ) , a,beC(A). (3.4) 



Recall that the Maurer-Cartan two-cocycle on a Lie algebra is a bilinear C-valued func- 
tion satisfying two conditions: 



Classical R-nmtrix theory of dispersionless systems: II 



5 



(i) it is skew-symmetric 

u 2 (a, b) = —0J 2 (b, a); (3.5) 

(ii) it satisfies the Jacobi identity 

u 2 ([a, b), c) + tu 2 ([c, a], b) + 0J 2 ([b, c], a) = 0. (3.6) 

Hence, the Lie product ( |3~3"D is well defined on C(A). The scalar product on C(A) is 
given by 

((a, a), {b, /3)) m := Tr(a • 6) + a ■ /3. (3.7) 
The Poisson bracket {-, •} on the functionals T>(C(A)) we define as 

{H,F}(L) :={(L,l),[(dF,l),(dH,l)]) m 

= ((L,l),([dF,dH),u; 2 (dF,dH)))c {A) , (3.8) 

for all (L, 1) e C(A*) = C(A). Then from Q we get the following form 

{H, F}(L) = (L, [dF, dH]) c(A) + u 2 (dF, dH), (3.9) 

which can be considered as a centrally extended Lie- Poisson bracket. 

Let us repeat the i?-matrix approach for the current Lie algebra C(A) with a natural 
Lie-Poisson bracket (|3i 



Lemma 3.1 Casimir functionals C G T>(C(A)) of a Lie-Poisson bracket ( |3.9|) satisfy 
the so-called Novikov-Lax equation 

[dC,L] + (dC) y = 0, (3.10) 

for allLeC(A*) ^C(A). 

Proof. For every H, F e T>(C(A)) and L G C(A) 

{H, F}(L) = (L, [dF, dH}) c{ A) + uJ 2 (dF } dH) 

= (dF, [dH, L}) c( a) + (dF, (dH) y }) c(A) = (dF, [dH, L] + (dH) y }) c(A) , 

hence for Casimir functionals C G T>(C(A)) 

{C, F}(L) = ^ [dC, L] + (dC) y = 0. □ 

The .R-structure R G End(C(A)) is defined as follows: 

[(a, a), (b,(3)]^ : .= ([a,b] R ,u*(a,b)), (3.11) 

where oo 2 (a, b) := u 2 (Ra,b) + uj 2 (a,Rb). Then, the new linear Lie-Poisson bracket has 
the following form 

{H, F}i(L) = ((L, 1), [(dF, 1), (dH, l)y ?(A) 

= (L, [dF, dH) R ) ei A) + u*(dF, dH). (3.12) 



Lemma 3.2 The following Poisson operator is related to the linear Poisson bracket 
6(L) : dH ^ -ad L RdH - R*ad L dH + (RdH) y + R*(dH) y . (3.13) 
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Proof. For every H , F G V(C(A)) and L G C(A) 

{H, F} l (L) = (L, [dF, dH] R ) c{A) + w*(dF, dH) 

= (RdF, [dH, L]) C(A) + {dF, [RdH, L]) C{A) + {RdF, (dH) y ) c(A ) 

+ (dF, (RdH)y)c(A) 
= (dF, —[L, RdH] - R*[L, dH] + (RdH) y + R*(dH) y ) c(A ) 
= (dF,9(L)dH) c{A) . □ 



Theorem 3.3 The Casimir junctionals C% G T>(C(A)) of the Poisson bracket (3.8) 
on C(A*) = C(A) are in involution with respect to the linear Poisson bracket ( S.li ). 
Moreover, Casimir junctionals Ci satisfy the following hierarchy of evolution equations: 



L ti = 6(L)dCi = [Rdd, L) + (Rdd) y , i G Z. (3.14) 

Proof. Let Cj and Cj G T>(C(A)) are Casimir functional, then 

{C t , Cj}i(L) = (L, [dCj, dCiUhw + ^(dCj, dQ) 
= (RdC 3 , [dC t , L) + (dC t ) y ) c(A) + (RdQ, [L, dC 3 ] - (dC^A) = 0. 

The proof of the second part of the theorem is obvious. □ 

In a special case of Poisson algebras, which are considered in the paper, the bracket 
( |3.12| ) is nothing else but a centrally extended Lie- Poisson bracket ( |2.9| ). For higher 
order Poisson brackets (|2.6| ) we failed to prove the Poisson property (Jacoby identity) 
after central extension. 



4 Poisson algebras of formal Laurent series 

Let A be an algebra of Laurent series with respect to p 

A=\l = Y j ^(x)A, (4.1) 

I i&1 J 

where the coefficients Ui(x) are smooth functions. It is obviously commutative and 
associative algebra under multiplication. The Lie-bracket can be introduced in infinitely 
many ways as 

dh x dh 2 dL x dL 2 

[Ll > L2] = P { ^p-^x~ ~ := L2} - r Z ' (4 ' 2) 

as adi = P r (^-§^ — %~§~) is a derivation of the multiplication, so A r := (A, {■, -} r ) are 
Poisson algebras. An appropriate symmetric product on A r is given by a trace form 
(a, b) A := tr(afo): 

trL = / res r Ldx, res r L = u r -i(x), (4.3) 
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which is ad- invariant. In expression ( |4.3| ) the integration denotes the equivalence class 
of differential expressions modulo total derivatives. For a given functional F(L) = 
J n f(u)dx, we define its gradient as 

where Sf/Sui is a variational derivative. 

We construct the simplest i?-matrix, through a decomposition of A into a direct sum 
of Lie subalgebras. For a fixed r let 



A^- r+k = P^- r+k A = < L = ^ u i( x )p l f 

A<- r +k = P<-r+kA = < L = \ 

I i<-r+k ) 



(4.5) 



where P are appropriate projections. As we presented in 0, A^ r+k , A^ r+ k are Lie 
subalgebras in the following cases: 

1. k = 0, r = 0, 

2. k = 1,2, r G Z, 

which one can see through a simple inspection. Then, the i?-matrix is given by the 
projections 

R = ^{P^-T+k — P<-r+k) — P^-r+k ~ 2 = 2 ~ ^ > <- r+k ' (^-6) 

To find i?* one has to find P^_ r+k and P<_ r+k given by the orthogonality relations 

(P^-r+ki P<-r+k) = (P<_ r+fc , P^-r+k) = 0. (4.7) 

So, we have 

P^-r+k = P<2r-k, P<-r+k = P^2r-k, (4.8) 

and then 

1 * * 1 1 

P* = 2~(P^-r+k - P<-r+k) = 2~ P^2r-k = P<2r-k ~ ~- (4.9) 

5 Centrally extended Poisson algebras of Laurent 
series 

Let A be an algebra of Laurent series with respect to p 

A=lL = Y,u l {x,y)A, (5.1) 

I i& J 
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where the coefficients Ui(x,y) are smooth functions of two variables x and y. As in 
(l+l)-dimensional case p was a conjugate coordinate related to x, let us now introduce 
q as a conjugate coordinate related to y. Then, introducing the extended Lie-bracket 
j|) in the form 

{Li L 2 } ■= p r ( — — - + ^1^2 _ 0LidL2 

dp dx dx dp dq dy dy dq ' 



and the extended Lax element £ = L — q, L £ A. The Lax-Novikov equation (|3.10| ) 
takes the form 

{dC,C}r = 0, (5.3) 

and the hierarchy of evolution equations ( |3.14| ) for Casimir functionals C(L) with R- 
matrix given by (|4.6j ) has the form of two equivalent representations 



C ti = {(dCi)^_ r+fe ,£} r = -{(dC i )<_ r . +jfc ,£} r , ieZ, (5.4) 

which are Lax hierarchies. 

To construct dispersionless (2+l)-dimensional integrable equations, at first we have 
to solve equation (|5.3|), which can be done by putting 



dC{ = CLjP^ , i ^ —r + k, (5.5) 



or by 



dCi = cijP' , i < —t + k, (5.6) 



where the function parameters dj are obtained from ( |5.3|) successively via the recurrent 
procedure. Notice that although the solutions ( |5.5| ) or ( |5.6| ) are in the form of infinite 
series, in fact we need only their finite parts (c?Cj)^_ r+ fc or (dCj) < _ r+ fe. Hence, for a 



given C in principle we can construct two different hierarchies of Lax equations (|5.4j ). 

We have to explain what type of Lax operator can be used in ( |5.4| ) to obtain a con- 
sistent operator evolution equivalent to some nonlinear integrable equation. Obviously, 
we are interested in extracting closed systems for a finite number of fields. Hence, we 
start with looking for Lax operators £ in the general form 

C = u N p N + u N ^p N - 1 + ... + u_ m+1 p- m+1 + u^ m p- m - q (5.7) 

of N-th order, parametrized by finite number of fields Ui. To obtain a consistent Lax 
equation, the Lax operator (|5.7|) has to form proper submanifold of the full Poisson 
algebra under consideration, i.e. the left and right-hand sides of expression ( |5.4|) have 
to lie inside of this submanifold. 

Observing ( |5.4|) with some (<iC) < _ r . + fc = a^ r+ k-\P~ r+k ~ l + a_ r+k ^ 2 p~ r+k ~ 2 + ■■■ one 
immediately obtains the highest order of the right-hand side of Lax equation as 

C t = {u N ) t p N + {u N ^) t p N - l + ... 

= -{(rfC)<_ r+ fc,M A rp iV + lower},. - d y (dC) < - r+k 
= (-((-r +k- l)a_ r . +fc _i(u A r) :!; - N (a„ r+k „i) x u N )p N+k ~ 2 + lower) 
+ (-(a-r+fe-i)^ r+fc_1 + lower) , (5.8) 
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where lower represents lower orders. Observing ( |5.4j ) with some (dC)^_ r+k = ... + 
a- r+k+ ip~ r+k+1 + a^ r+k p~ r+k one immediately obtains the lowest order of the right- 



hand side of Lax equation ( |5.4| ) as 

C t = ... + (u. m+1 ) t p- m+1 + (u_ m ) t p- m 

= {(dC)^- r+k , higher + U- m p~ m } r + d y (dC)^- r+k 

= (higher + ((-r + k)a_ r+k {u^ m ) x - (-m)(a_ r+fe ) :r M_ m )p" m+fc ' 1 ) 

+ (higher + (a_ r+fc ) v p~ r+fc ) , (5.9) 



where higher represents higher orders. Simple consideration of ( |5.8| ) and ( |5.9| ) with 
condition N ^ — m leads to the admissible Lax polynomials with a finite number of 
field coordinates, which form proper submanifolds of Poisson subalgebras. They are 
given in the form 

k = 0, r = : 

C = c N p N + c N ^p N ~ l + u N - 2 p N - 2 + ... + u -q for N ^ 1, (5.10) 
k = l,r e Z: 

£ = c N p N + UN _ 1 p N ~ 1 + „, + u _ m p- m - q f or N^l-r^-m, (5.11) 
£ = u^ r p~ r + M__ r _ip~ r_1 + ... + u^ m p~ m — q for — r ^ — m, (5-12) 
k — 2, r e Z : 

C = u N p N + ... + u 1 ^ m p 1 - m + c. m p- m -q for N ^ 1 - r ^ -m, (5.13) 
C = u N p N + ... + u 3 - r p 3 - r + u 2 - r p 2 ~ r - q for iV^2-r, (5.14) 

where the Ui are dynamical fields and cn, cn-i, c_ m are arbitrary time independent 
functions of x and Lax operators for k = 0, 1,2: ( |5.iq ),( |5TLT| ),( |5.13"|) form a proper 



submanifold in (l+l)-dimension 0, hence the Lax dynamics induced by them can be 
reduced onto (l+l)-dimensional space. Lax operators for k — 1: ( [5.1 2| ) and k = 2: 
(|5.14j) do not form s proper submanifold in (l+l)-dimension, hence the Lax dynamic 
induced by them is purely (2+l)-dimensional effect, and they cannot be reduced onto 
(l + l)-dimensional space. 

Hence, we know the restricted Lax operators £ we can now investigate the form of 
gradients of Casimir functionals dCi given by ( |5.5|) or by ( |5.6| ) which satisfy equation 
(|5.3| ), as well as we can investigate some further simplest admissible reductions of Lax 
operators. 



The case of k — 0. Let us consider Lax operators of the form (|5.10| ). Then observing 



|) with some dCi = dip 1 + a^ip 1 1 + lower one immediately obtains the conditions 
for the highest terms of dCi, since 

{a i p i + ai^p^ 1 + a^P^ 2 + lower, £} = -N(a i ) x c N p i+N ~ 1 

- (Nc N ( ai ) x + (N — ^CN^ai^^p 1 ^ 2 + lower = 0. (5.15) 



Therefore (aj) x . = (cij-i)a; = 0, idi(uN-2)x — NcN(a>i-2)x = and so on, hence (|5.5| ) has 
the following form 

dd = dip 1 + cti-ip l ~ l + —^—un^P 1 ' 2 + ai-sP l ~ 3 + lower, i ^ 0, (5.16) 

Nc N 
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where a,, are arbitrary x independent functions. Observing ( |5.3| ) with dCi ■ 
higher + di + ip l+1 + dip 1 one obtains the conditions for the lowest terms of dCi, since 



{higher + a i+ ip l+ + a t p l , C} = 

higher + ((i + l)a i+1 (u ) x + ia i {u 1 ) x - (a^Ui + {a i ) y )p t + iai(u ) x p 



i-1 



0. 



(5.17) 

Accordingly aj = and <Zj_i = <2j_ 2 = ... = since dj depends linearly on a J+1 , ctj +2 , ■ a%- 
Hence for k = there is only one Lax hierarchy for the dCi of the form ( 5.16 ). There 
are no any obvious further reductions of C. 



The case of k = 1 . For Lax operators of the form (|5.11|) by observing ( |5.3|) , dCj given 



by ( f5.5|) or (|5.6|) have the following forms 



dCj = «jj9* H — mat-iJ5* _1 + a^ 2 p l ~ 2 + lower, 

Nc N 

dCi = higher + a i+2 p l+2 + a i+ ip l+1 + aj(u_ m ) _ 



~p' 



i — r + 1, 
i < — r + 1, 



(5.18) 
(5.19) 



where «j is an arbitrary x-independent function. For Lax operators of the form ( 5.12 ) 
by observing (|5.3|), dCi given by ( |5T5] ) or ( |5.(j| ) have the following forms 



cZCj = /5y - 9^ 1 (ifi i (u- r ) x + r((3i) x u^ r )p l 1 + a^ 2 p l 2 + lower, i > -r + 1, 

(5.20) 

dCj = /li^/ier + a i+2 p i+2 + a i+ ip i+1 + ai(M_ m ) - -^, i < -r + 1, (5.21) 



where aj and /3, are arbitrary x- and y-independent functions, respectively. 

We remark that, if — m < 1 — r in C, there is a further admissible reduction of the 
equations (|5.4j), given by w_ m = 0, since such reduced Lax polynomials still are of the 
form (|5.11 ) or ( 5.12|) . We have to look for the form of gradients of Casimir functionals 
after such a reduction. It is easy to see that by this reduction u_ m = 0, the gradients of 
Casimir functionals (|5.18|) and ( |5.20|) preserve the order of the highest terms, and the 
form. For gradients of Casimir functionals Q5.19 ) and ( 5.21Q by this reduction the lowest 
order disappear, and as all other terms depend linearly on it, such gradients reduce to 
zero, except the one case (<iC i ) < _ r+ i = (L)<_ r+1 which produces equation £ t . = —C y . 
We can continue the reductions by putting «i_ m = 0, if the reduced £ are still of the 
form ( 5.11 ) or Q5.12 ) and so on. Therefore, the reductions are proper in general only for 
the gradients of Casimir functionals in the form ( |5 . 1 8 ) and ( |5.20| ). 



The case of k = 2. For Lax operators of the form (|5.13|) by observing (|5.3[) , dCi given 



by (p-5|) or (5.6) have the following form 



dCi = ai{u^) N p % + cii-ip 1 + di-2P % + lower 



dCi = higher + cii +2 p 



i+2 



lOti 



U X - m p lJrl + OLip 1 



mc. 



i ^ -r + 2, 
i < -r + 2, 



(5.22) 
(5.23) 
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where oci is arbitrary x-independent function. For Lax operators of the form ( |5.14| ) by 
observing ( |5.3|) , dCi given by ( |5.5|) or (|5.6|) have the following form 



oci(u N ) * N p l + di-ip 1 1 + di-2V % 2 + lower 



i ^ -r + 2, 

higher - ar 1 (i/3 i (ii2-r)x - (2 - r)(ft),« 2 -r)p !+1 + Ap\ 



z < — r 



(5.24) 
-2, 

(5.25) 



where a* and /3, are arbitrary x- and y-independent functions, respectively. 

If TV > 1 — r in C, there is a further admissible reduction of equations (|5.4|) , given 
by un = since such reduced Lax polynomials still are of the form Q5.13|) or (|5.14j ). By 
analogous considerations as for k = 1, these reductions are proper in general only for 
the gradients of Casimir functionals in the form ( |5.23[) and (|5.25f) . 

The different schemes are interrelated as it is explained in the following theorem. 



Theorem 5.1 Under the transformation 
x, y = -y, p = p , q = q, t 



x 



t 



(5.26) 



the Lax hierarchy defined by k = 1, r and C transforms into the Lax hierarchy defined 
by k = 2, r' = 2 — r and C , i.e. 



k — 1, r, C 



k = 2, r' = 2-r, C. 



(5.27) 



Proof. It is readily seen that the Lax operators for k = 1 and r of the forms ( |5.11| ), 
Q5.12 ) transform into the well restricted Lax operators for k = 2 and r' = 2 — r of the 
forms ( 5.13|) , ( |5.14j ) respectively. Let's observe that 



_ pr ^AdB_ _ dAdB_ + dAdB_ _ dAdB_ 
' r dp dx dx dp dq dy dy dq 



i-r+2 

= —p 

and 

(dcy >s = (dC)^ s 

Hence, we have 
Ct 



dA'dB' dA'dB' 



dp' dx' dx' dp' dq' dy' dy' dq' 



{(dc)^ r+1 , c} r = -{{dc)^_ r+u cy r , 

— {(dC )^ r _ 1( £ } r , = —{(rfC / ) < _ r / +2 , £'}',, = C t ,. 



□ 



Therefore, some dispersionless systems can be reconstructed from different Poisson 
algebras. Moreover, we remark that the gradients of Casimir functionals for k = 1 and 
k = 2 transform by p^ 1 = p' reciprocally at slant, i.e. flgl8D <-» (fT23D , d5l9D <-> (pT24D and 
(|5T20|) ^(|5T25| ) , (|5T2T|) ^(|5^ ) . 

Two equivalent representations of Poisson structure coming from the linear Poisson 
tensor ( |3.13| ) with the i?-matrix given by (|4.6|) are 



6{C)dH = {(dH)^ r+k ,C} r - ({dH,C} r )^ 2r _ k 

= — {(rfiJ) < _ r+fc , £} r + ({dH, £} r )<2r-A;- 



(5.28) 
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It turns out that the first representation yields a direct access to the lowest polynomial 
order of 6dH, whereas the second representation yields the information about the highest 
orders present. There are two options. The best situation is when a given Lax operator 
forms a proper submanifold of the full Poisson algebra, i.e. the image of the Poisson 
operator 9 lie in the space tangent to this submanifold for each element. If this is not 
the case, the Dime reduction can be invoked for restriction of a given Poisson tensor to 
a suitable submanifold. 



The case of k = 0. Let us first consider the simplest admissible Lax polynomial ( 5.1C ) 
of the form 



£ = p N + u N _ 2 p N 2 + ••• + uip + u 



q- 



(5.29) 



This is the well-known dispersionless Gelfand-Dickey case. Then, the gradient of the 
functional H(C) is given in the form 



8_H_ 

Jc 



5H 

5w 



-P 



5H 

5ui 



P 



5H 

8u N - 



-P 



A-N 



(5.30) 



By inserting ( |5.29| ) into ( |5.28| ) it becomes clear from the first representation of the 
linear tensor that lowest order of 6dH is at least zero, from the second representation 
it is evident that the highest differential order will be at most N — 2. Hence, 9dH is 
tangent to the submanifold formed by the Lax operator of the form ( |5.29| ). As a result, 
these Lax operators form a proper submanifold of full Poisson algebra, and the Poisson 
tensor, since (§^) >0 = 0, is given by 



4 



5H\ 



(5.31) 



o/ >o 



The case of k 
in the form 



1. Let us first consider the simplest admissible Lax operator ( |5.11| ) 

(5.32) 



C = p N + u N ^p N ~ l + ... + u^ m p x ' m + u. m p~ m - q 



Then gradient of functional H(C) is given in the form 



6H_ 



_5H_ 



y+m-1 + 



SH 



_pr+m-2 _|_ 



m+1 



SH 



J--N 



-p' 



(5.33) 



Inserting ( ^32|) into ( ggp one immediately obtains the highest and lowest order of 8dH 



as 



9dH = ((-.y^ 1 + lower) + ((..> 2r ~ 2 + lower) 
= (higher + (...)p~ m ) + (higher + (...)p 2r_1 ) 



(5.34) 



where lower (higher) represents lower (higher) orders. Hence, Lax operators of the form 
(|5.32|) form a proper submanifold for N ^ 2r — 1 ^ —m, as then 6dH is tangent to this 
submanifold. So the linear Poisson map is 



•(g) 




£ + 



c 



SH_ 



(5.35) 



r/ >2r-l 
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Otherwise a Dirac reduction is required. 

Analogously, for Lax operators (|5.12| ) in the form 

£ = u^ r p' r + U- r -ip~ r ~ X + ... + u^p 1 "™ + u. m p~ m - q (5.36) 

we have 

OdH = ((...)p" r + lower) + ((...)p 2r ~ 2 + lower) 

= (higher + (...)p~ m ) + (higher + {...)p 2r ' 1 ) . (5.37) 

Hence, this operator forms a proper submanifold for r ^ and 2r — 1 ^ —m. The 
Poisson tensor is given by ( |5.35| ). In other cases a Dirac reduction is required. The 
simplest case is r = 1 with one-field reduction. Let 

C = u + C = u + u^ip^ 1 + u- 2 p~ 2 + ••• + u X - m p X ~ m + u- m p- m - q. (5.38) 

The Dirac reduction with the constraint u = leads to the Poisson map in the form 

which is generally nonlocal. 



The case of k = 2. Let us consider Lax polynomials ( |5.13|) in the form 

C = u N p N + u N ^ lP N - 1 + ... + u^mp 1 -" 1 + p~ m - q. (5.40) 
Then gradient of functional H(C) is given in the form 

^ = ^P r+m ~ 2 + - + 7—— — p r ~ N + ^Lf- N ~\ (5.41) 

6L 0Ui- m OUn-1 OUn 



Then by analogous considerations as for k = 1 or by Theorem |5.1| , C given by ( |5.40| ) 



form a proper submanifold for A^ ^ 2r — 3 ^ —m. The Poisson tensor has the form 

Otherwise a Dirac reduction is required. 

Analogously, Lax operators ( |5.14| ) in the form 

C = U N p N + UN-^' 1 + ... + M 3 -rP 3 ~ r + U 2 -rP 2 ~ r ~ q, (5.43) 

form a proper submanifold for r ^ 2 and A" ^ 2r — 3. Then, the Poisson tensor has 
the form ( |5.42j ). Otherwise a Dirac reduction is required. The simplest case is for r = 1 



with one-field reduction. Let 

C = £ + u = u N p N + u N ^ip N ^ + ... + u 2 p 2 + uip + u — q. (5.44) 
The Dirac reduction with the constraint u = leads to the Poisson map in the form 



^(iM{*£}l-W*£H- 



(5.45) 
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which is generally nonlocal. 

Hence, we know the Poisson structure for (2+l)-dispersionless systems constructed 
from Poisson algebras, and since we are interested in Hamiltonian systems, we shall now 
consider the problem of their construction. The conserved quantities Hi are described 
by the Hamiltonian equations 

C ti = OdHiiC). (5.46) 

First we have to find cosymmetries (one-forms) dHi which are gradients of Hamiltonians. 
Because we are using gradients of Casimir functionals dCi to generate equations (|5.4p, 
our dHi are given by projections of dCi on subspaces spanned by dHi in the form 
( [5.30|) , Q5.33 ) and Q5.41 ) for k — 0, 1 and 2, respectively. Then, we can apply the Lemma 



6| and hence Hamiltonians are defined as follows 

1 r r rl 



Hi(C) = [ Tr(dHi(XC)C)dX = [[ [ res r (dHi(\£)£)d\dxdy. (5.47) 
Jo J Jnxs 1 Jo 

For Lax operator £ = Yli=i u iV l ~ Q the gradients from ( |5.47|) are given by 



- fih 

dH t (XC) = V 77r ^(An 1 , Au 2 , \u n )p r ~ l ~\ (5.48) 
Hence, by using definition of the residuum ( |4.4| ) we get that 

res r (dHi(\C)C) = ^ Ui J(\^^^ Ul, ' Ku ' 2 ' ""' ' KUn ^ (5.49) 



i=l 



Contrary to the (1+1) -dimensional case in the (2+1) case, the functional densities 
contain terms with x, y derivatives as well as nonlocal terms. Nevertheless, all these 
additional terms appear in a special form, namely they are expressed through the pseudo- 
differential operators of the form U k ,U~ k where 

V:=d~ x d y , V 1 : i) u -i) r . (5.50) 

Thus, the additional to Q5.47 ) useful relation of calculation of variations containing V, 



derived from Q2.2|) , is the following one 

6 H f{u)V k g{u) dxdy = ^^V k g{u) + ^^-V k f(u). (5.51) 



JJ nxS i du du 



6 A list of some (2+l)-dimensional dispersionless 
systems 

In this section we will display a list of the simplest nonlinear dispersionless (2+1)- 
dimensional integrable systems. Calculating the gradients dC n (n-highest order) given 
by ( |5.5| ) we consider the Lax hierarchy 

C tn = {(dC n )^_ r+k ,C} r , neZ. (6.1) 
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The second hierarchy for dC n given by (|5.6|) can be obtained by the transformation from 
Theorem |5.1| , which we leave for the interested reader. We present the Hamiltonian 
structure for particular choices of r. For k = and k = 1 the choice n = 1 — r will 
always lead to the dynamics (uijt^ = (1 — r)(ui) x for the fields Ui in L, so that we may 
identify ti__ r = jz^x in this cases. For (dC n ) j>_ r +fc = L the equations become trivial, 
and then C t • , = £v For each choice of k = 0, 1 or 2 and N we will exhibit the first 

L n-\-r — k y ' 

nontrivial of the nonlinear Lax equations (|6.1|) associated with a chosen operator C. 



The case of k — 0. 

Example 6.1 The dispersionless Kadomptsev-Petviashvili: k = 0, r = 0, iV = 2. 

The dispersionless Kadomptsev-Petviashvili (dKP) equation is a (2+l)-dimensional 
extension of the dispersionless KdV equation. The Lax operator for the (2+l)-dimensional 
dKP hierarchy has the form 

£ = p 2 + u - q. (6.2) 

Then we derive for {dC^)^o = p 3 + \up + \T>u 

3 3 

u t3 = -uu x + -Vu y = 9dH, (6.3) 
where we get the Poisson tensor and the Hamiltonian 

9 = 2d x , H = - I (w 3 + -uV 2 u) dxdy. (6.4) 

Example 6.2 The (2+1) Boussinesq hierarchy: k = 0, r = 0, iV = 3. 
The Lax operator is given by 

C = p 3 + up + v — q. (6.5) 

We derive for (dC 2 )^o = p 2 + \u 



3 

" \ ( 2v x 

t2 ~ \ ~\ uu * + l u v 



9dH. (6.6) 



Eliminating the field v from this equation we can derive the (2+l)-dimensional 'disper- 
sionless' Boussinesq equation 

4 2 

u tt = -^u xy - -{u 2 ) xx . (6.7) 
The respective Poisson tensor and Hamiltonian are given in the following form 

3d, \ „ Iff , 1 3 , 2,1 



6 = 1 30, o* j • H = 3 JLJ-r + " + r Vu) dxdy - (6 - 8) 
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Example 6.3 The case: k = 0, r = 0, iV = 4. 
The Lax operator is 

C = p 4 + up 2 + vp + w — q, (6.9) 
then for (dC 2 )^o = p 2 + \u we have 

2v x 

2w x -uu x | = 6dH, (6.10) 



V 













where 



o o 4a 

4^ 0' j. (6.1.1.) 
4<9 cLw + ud 



H = - [ {~\u 2 v + 2vw + -uVu) dxdy. (6.12) 
The case of /c = 1. 

Example 6.4 Three field hierarchy: fc = l,r£Z\2. 

The Lax operator has the form (|5.11|) with N = 2 — r, m = r + l 



£ = p 2 ~ r + up 1 ' 7 ' + vp- r + wp' 1 " 1 - q. (6.13) 
Then for (dC2- r )^-r+i — P 2 ~ r + up l ~ r we have 

u y + (2 - r)u x \ 
ru x t> + (1 — r)w, + (2 — r)w x . (6.14) 
(1 + r)u x w + (1 — r)uw x J 

This Lax operator forms a proper submanifold as regards the condition N ^ 2r— 1 ^ — m 
only for r = 0, 1, otherwise a Dirac reduction is required. Then for r = 

My + 2V X 

uv x + 2w x I = fldF, (6.15) 







• 






*2-r 



V 













where 



2a„ 



29, u^-d, I , (6.16) 



2d x d x u-d y 



H = — I [ (lGvw - 2u 2 Vv + 8uVw + -u 2 Vu 2 
lGJJtix&\ 1 



+4v£>i; - mD^ 2 + AuV 2 v + uV 6 u dxdy. (6.17) 
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For r = 1 we have 



V 




( Uy + V X \ 








u x v + w x 


= 6dH, 


w 




V 2u x w J 





(6.18) 



where 



d y d x v , ! 



| ^ d x w + wd x , H = - (u 2 + 2v)dxdy. (6.19) 

2^ / 2 JJnx& 

Example 6.5 Dispersionless (2+1) Toda: k = l,r G Z\ {2}. 



The first admissible reduction u; = of (|6.13|) leads to the two field Lax operator 

C = p 2 ~ r + up 1 "*' + vp~ r - q. (6.20) 

This Lax operator forms a proper submanifold only for r = 1, otherwise a Dirac reduc- 
tion is required. Hence, for r = 1 by reduction w = ( 6.18| ) we get the (2+l)-dimensional 
dispersionless Toda equation 



u \ ( u y + v x 
v I , \ u T v 



6dH, (6.21) 



where 





( 4 a o)- H -lILy +2v)dxdv ' < e22 > 



known till now in a few non-Hamiltonian representations [l], |2|, l7fl . Changing the inde- 
pendent coordinate t' = t — y and eliminating w-field one gets 

(lnv)w = v xx or (j) W = (e 0a % (6.23) 

where d> v = In v. For r = we have 



12 



u \ I u y + 2v x 



(6.24) 



but we lose the Hamiltonian structure since the Poisson tensor ( |6.16| ) is not reducible 
with the constraint w = 0. Hence, the Lax operator (|6.20|) for r = generates equations 
which are non-Hamiltonian. 

The next admissible reduction w = v = of ( |6.14| ) leads to trivial equation L t2 _ r = 
L y since (dC2- r )^- r +i = L- 

Example 6.6 One field hierarchy: k = l,r G Z \ {2}. 
The Lax operator is given in the form 



£ = p 2 ~ r + (2 - r) V~ r - q- 



(6.25) 
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Then one finds for {dC 3 . r )^_ r+1 = p 3 ~ r + (3 - r)up 2 - r + (^Vu + |(3 - r)u 2 ) p l ~ r a 
whole family of (2+l)-dimensional dispersionless one-field systems 

u t3 _ r = ~{3-r)(l-r)u 2 u x + T ^~^ uu y + t|^2 ^%+ - ~ ^ u x Vu, (6.26) 



derived for the first time in fllql , including the modified dKP as a special case of r = 0. 

This Lax operator forms a proper submanifold only for r = 1, in other Dirac 
reduction is required. For r = 1 we get 

u t2 = 2uu y + 2Vu y = 6dH, (6.27) 

where 

6 = d y , H=[[ + uVu) dxdy. (6.28) 

For r = we get 

3 3 3 

u h = --u 2 u x + -T>u y + -u x Vu, (6.29) 



and by Dirac reduction of (|6.16 ) with the constraint w = v = we get the formal 
Poisson tensor 

6 red = 8d x (d y - 2ud x y 1 d x (d y - 2d x u)- l d x , (6.30) 
and the related sympletic tensor 

J = {e^y 1 = ±(2> - 2u)d7 x \V - 2tt), (6.31) 
such that Ju t3 = dH, where 

H=?-[[ (-\u & + uVu A + \u 2 V 2 u 2 - u 2 (Vu) 2 
32 J Jnxs 1 V 3 2 

+ ^(Vu) 3 - uV 3 u 2 + ^uV 4 t?J dxdy. (6.32) 

Example 6.7 Three field hierarchy: k = 1, r £ Z. 

This case does not exist in (1 + 1) -dimension. The Lax operator has the form ( 5.12 ) 
with m = r + 2 

£ = -u£r r + vpT r ~ x + wp~ r ~ 2 — g. (6.33) 
Then for (dC2- r )^-r+i = p 2_r + (r — 2)'D~ 1 u]9 1 ~ r we have 



u 

V 

w 



i-2-r 



ruD l u x + (1 — r)u x T> l u — v x 
(r-2)( {l+r)vV- l u x + {l-r)v x V- 1 u-w x |. (6.34) 
(2 + r)wV- l u x + (1 - r)w x V- x u 
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This Lax operator forms a proper submanifold as regards the condition 2r — 1 ^ — m 
only for r = 0, otherwise a Dirac reduction is required. Then for r = 







V 




w 





u x T>~~ x u — v x 
T>~ x u x + v x T)~ x u — u> 
2wT>~ 1 u x + w x T>~ x u 



8dH, (6.35) 



where 






I , (6.36) 



c^tu + tyc^ 



and 

H 



For r = 1 we have 



^uP-^-vZrWi^Xr 1 ^) 2 ) dsdy. (6.37) 







• 













■uZ> l u x — v x 

2vV~ 1 u x - w x | = 0d£f. (6.38) 
?>vfD~ x u x 

We derive the Poisson tensor from ( |5.39| ), then 

-uD~ 1 d x u + d x v + vd x —2uD~ x d x v + 2d x w + wd x —3uV~ 1 d x w 

- "' 1 ~2vV~ x d x u + d x w + 2wd x -AvV~ l d x v -QvV^d^ . . 



-3wV 1 d x u —6wV 1 d x v —9wV 1 d x w 



(6.39) 



and 



H = u dxdy. (6.40) 

J JQxS 1 

Example 6.8 Two field hierarchy: k = l,r 6 Z. 

The first admissible reduction w = of (|6.33|) leads to the two field Lax operator 

C = up~ r + vp- r ~ l - q. (6.41) 

This Lax operator forms a proper submanifold only for r = 0, otherwise a Dirac reduc- 
tion is required. Hence, for r = by reduction w = of (|6.35|) we get 



i; ^ " y vT>^ l u x + v x T>^ L u 



where 
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H= - It (-vV- l v + v (V^u) 2 ) dxdy. (6.44) 
For r = 1 we have 

:) ti =-("i;x;0=^ (M5) 

We derive the Poisson tensor from ( |6.39| ) with the constraint w — 0, then 

nred _ ( -uV~ l d x U + X V + V& x -1uV~ X d x V \ 

U ~\ -1vV- x d x u -AvV-^v '' ' ' 



and 



H = u dxdy. (6.47) 

Example 6.9 One field hierarchy: k = l,r G Z. 

The second admissible reduction u> = u = of (|6.33|) leads to the one field Lax 
operator 

C = up- r - q. (6.48) 

This Lax operator does not form a proper submanifold as the condition 2r — 1 ^ — m 
is violated, hence a Dirac reduction is required. For r = by reduction v = w = of 
(|6.35|) we get 

u t2 = -2u x V-\ (6.49) 

but we lose the Hamiltonian structure as the Poisson tensor Q6.43 ) is not Dirac reducible 
with constraint v = w = 0. For r = 1 we have 

u h = -uD~ l u x = 9dH. (6.50) 

We derive the Poisson tensor from (|6.46|) with the constraint w = 0, then 

9 red = -uD- x d x u, (6.51) 

and 



H = u dxdy. (6.52) 

J JQxS 1 

The case of k = 2. 

Example 6.10 One field hierarchy: k = 2, r G Z \ {2}. 
The simplest admissible Lax operator is given by 

C = u 2 ~ T p 2 - r - q. (6.53) 

This case does not exist in (l+l)-dimension. In this case we have to consider separately 
two cases: r ^ 1 and r = 1. Then, one finds again a whole family of (2+l)-dimensional 
dispersionless one-field systems |16|] including a dispersionless (2+l)-dimensional Harry 
Dym equation as a special case of r = 0: 



Classical R-nmtrix theory of dispersionless systems: II 



21 



r^l : for (rfC 3 _ r .)^_ r+2 = u 3 r p 3 r + — fy^r^yW 2 r 'Du r 1 p 2 r one finds 



3 3 7* 

r = 1 : for (dC 2 ) >1 = u 2 p 2 + 2uD (In it) p one finds 

w t2 = 2u y V In m + 2w£> (In u) y . (6.55) 

To get 9, we have to make a Dirac reduction as the conditions r ^ 2, iV ^ 2r — 3 
are violated. Poisson tensor for r = 1 is given by (|5.45 ), then we get for (|6.55|) the 
Hamiltonian structure, where 

9 red = uV- l d x u, H=ll (in uD 3 In m + — (Vmuf ) dxdy. (6.56) 

Example 6.11 Two /ieW hierarchy: k = 2, r e Z \ {3}. 
The Lax operator is given by 

L = V~ r + - q- (6.57) 

This case is nonreducible to (l+l)-dimension. Then, one finds for ( k dC2~ r )>- r +2 — 
u 3 ~ r p r 

u\ _ 2 - r u _± / (3 - r)n^ - (2 - r)^ \ ^ 



To get # we have to make a Dirac reduction as the conditions r ^ 2, iV ^ 2r — 3 are 
violated. The Poisson tensor for r = 1 is given by ( |5.45| ), then 



u \ v u I 2uv x — u x v 



9 red dH, (6.59) 



where 



1 2vV o, r u 2vV y ajv + o^u + mo t ' 



and 



H=f[ - inP^ + ^Pln^ dxdy. (6.61) 

Example 6.12 Two field hierarchy: k = 2, r e Z \ {2}. 
The Lax operator is given by 

£ = u 2 - r p 2 - r + vp l ~ r + p- r - q. (6.62) 

Then, for {dC 2 - r )^-r+2 = u 2 ~ r p 2 ~ r one finds 

w ^ _ ( u y - (1 - r)u^ + w a ^ , g 



D / V (2 — r)ru r u a 

/ T2 — r 
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This Lax operator forms a proper submanifold only for r = 1, otherwise a Dirac reduc- 
tion is required. Hence, for r = 1 we get 



u \ I u y + uv x 
v . \ u~ 



6dH, (6.64) 



where 



and 



em = ( ^ | , (6.65) 



H = jj ( u + -v 2 + vT>\nu + - lnw V 2 \nu J cfrrdy. (6.66) 

Example 6.13 Three field hierarchy: k = 2, r G Z. 
The Lax operator is given by 

£ = wp 2_r + fp 1_r + u;p~ r + p~ r_1 — g. (6.67) 

Then for (c?C2- r )^_ r +2 = up 2_r one finds 

u y - (1- r)u x v + (2 - r)uv x \ 

ru x .w + (2 — r)uw x . (6.68) 

(1+rK / 

This Lax operator forms a proper submanifold only for r = 1, otherwise a Dirac reduc- 
tion is required. Hence, for r = 1 we get 

Uy + UV X 

u x w + uw x ] = 6dH } (6.69) 
2u~ 







[V 






t2-r 







• 











where 



and 



m9 x 

fV = J 9^ -d y | , (6.70) 
2^ 



H = / I ( —v 2 + uw + v T> In u H — hiMT> 2 lnu) dxdy. (6-71) 
'Qxs 1 \2 2 
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